A Sufficient Condition for Partial Ensemble Controllability of Bilinear 
Schrodinger Equations with Bounded Coupling Terms 
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Abstract — This note presents a sufficient condition for par- 
tial approximate ensemble controllability of a set of bilinear 
conservative systems in an infinite dimensional Hilbert space. 
The proof relies on classical geometric and averaging control 
techniques applied on finite dimensional approximation of the 
infinite dimensional system. The results are illustrated with the 
planar rotation of a linear molecule. 

I. INTRODUCTION 

A. Control of quantum systems 

The state of a quantum system evolving in a Riemannian 
manifold is described by its wave function, a point rp 
in L 2 (fl,C). When the system is submitted to an electric 
field (e.g., a laser), the time evolution of the wave function 
is given, under the dipolar approximation and neglecting 
decoherence, by the Schrodinger bilinear equation: 



.dip 



■A + V(x))ip(x, t) + u(t)W(x)i>(x, t) (1) 



where A is the Laplace-Beltrami operator on f2, V and W are 
real potential accounting for the properties of the free system 
and the control field respectively, while the real function of 
the time u accounts for the intensity of the laser. 

In view of applications (for instance in NMR), it is 
important to know whether and how it is possible to chose 
a suitable control u : [0, T] — > R in order to steer ([TJ 
from a given initial state to a given target. This question 
has raised considerable interest in the community in the last 
decade. After the negative results of [1] and [2] excluding 
exact controllability on the natural domain of the operator 
— A + V when W is bounded, the first, and at this day the 
only one, description of the attainable set for an example 
of bilinear quantum system was obtained by Beauchard ([3], 
[4]). Further investigations of the approximate controllability 
of ((TJ were conducted using Lyapunov techniques ([5], [6], 
[7], [8], [9], [10]) and geometric techniques ([11], [12]). 

B. Ensemble controllability 

In many applications, a macroscopic device acts on a 
large number of identical microscopic quantum systems (for 
instance, a single laser acts on a small quantity of liquid 
containing many molecules). Usually, the external field acts 
differently on each of the small systems (depending for 
instance on the orientation of the molecule with respect to 
the exterior electric field). For the sake of simplicity, we will 



assume in this work that the action of the external field on 
the system a is proportional to a. Instead of one system of 
type ([TJ, one has to control a continuum: 



.dipa 



= (-A + V(x)+u(t)aW(x))^ a (x,t) (2) 



where the system labeled with a, a £ [0, 1], has wave 
function i/j a . Notice that, since the systems are physically 
identical, the free dynamics (when u = 0) is the same for 
every a. 

The simultaneous (or ensemble) control problem turns into 
the following question: let a continuum of initial conditions 
W>2)ae[o,i] and of targets (^*) oe[ o,i] be given. Does it exist 
a control u that steers the systems (|2]l from ip® to ?/>* for 
every a in [0, 1]? 

Because of its crucial importance for applications (disper- 
sion of parameters is always present in real world systems), 
the problem of ensemble controllability of quantum systems 
has been tackled by many authors, see for instance [13], 
[14], [15], [16], [17], [18] for theoretical results and [19] for 
numerical aspects. 

C. Framework and notations 

To take advantage of the powerful tools of linear operators, 
we will reformulate the problem (|2) in the following abstract 
setting. Let H be a separable Hilbert space, endowed with the 
(•, •) Hilbert product. We consider the continuum of control 
systems 



^ A I 



u(t)aBip, a £ [0, 1] 



(3) 



where the linear operators A, B can be completed in a 5- or 
6-uple that satisfies Assumption Q] or Assumption [2] 
Assumption 1: The 6-uple (A, B, Ai, <j)\, A2, 4>2) satisfies 

1) A is skew-adjoint with domain D(A); 

2) B is bounded and skew-symmetric; 

3) 4>\ and <f>2 are two eigenvectors of A of norm 1, 
associated with eigenvalues — iAi and — iA2; 

4) Ai < A 2 and (</> u Bfa) ^ 0; 

5) for every eigenvalues — i/j, and — i/x' of A, associated 
with eigenvectors v and v' , |Ai — A2I = \fJ, — 
implies {Ai,A 2 } = {n, //} or {Ai,A 2 } n {/i,//} = 
or (v, Bv') = 0; 



6) the essential spectrum of LA (if any) does not accumu- 
late in any of these four points: 2Ai — A2, Ai, A2 and 
2A 2 -Ai. 

Assumption 2: The 5-uple (A, B, U, A, <&) satisfies 

1) A is skew-adjoint with domain D(A); 

2) B is skew-symmetric; 

3) U is a subset of R containing at least and the points 
{1/n, n G N}; 

4) for every u in U, A + uB is skew-adjoint (with domain 
not necessarily equal to D(A)); 

5) $ = (<fij)j e -N is a Hilbert basis of H made of 
eigenvectors of A, all of which in the domain of B; 

6) A = (Aj)jgN is a sequence of real numbers such that, 
for every j in N, A<j)j = —iXj; 

7) Ai < A 2 and (faBfa) ^ 0; 

8) for all eigenvalues — ifi and — ipj of A, associated with 
eigenvectors v and u', |Ai — A 2 | = \fi — fi'\ implies 
{Ai,A 2 } = {m,m'I or {Ai,A 2 } n {/i,//} = or 
(u,Bi;'> = 0. 

If (A, B) satisfies Assumptions [TTT1 and [Tl2l (resp. (A, £,£/) 
satisfies Assumptions 121 II 12131 and \2\A\ . then for every to,t 
in R, for every u in L 1 (R, R) (resp. u : R — > U piecewise 
constant), there exists a unique family of unitary operators 
( T "Oae[o,i] such that, for every family (V£)ae[o,i] in H > 
for every a in [0, 1], t i-» Y"' t "-0° is the unique solution of 
(0 in the weak sense that satisfies T"'" ip® = ip^,. 

Let (A, B, Ai, <j>i, A2, <fo) satisfy Assumption Q] or 
(A, B,U, A,(p) satisfy Assumption |2] We define the 2- 
dimensional Hilbert space £ 2 = span(0i, 2 ), 7r 2 : H — » 
the orthogonal projection on £ 2 , A' 2 ) = 7r 2 A7r 2 and £?( 2 ) = 
7r 2 A7r 2 the compressions of A and £? on £ 2 , and AT^'", the 
propagator associated with the (infinite dimensional) system 
x' = A^x + uaB&'x. By abuse, we will still denote with 
X^" the restriction of X^ to £ 2 . 

D. Main result 

Proposition 1: Let (A, B, Ai, <j>\, A 2 , 2 ) satisfy Assump- 
tion Q] (resp. (^4, B, U, A, $) satisfy Assumption |2]) and let 
T : a G [0, 1] h-> T a G £/(£ 2 ) be a continuous curve of 
unitary operators on £ 2 that satisfies T° = Idx 2 . Then, for 
every e > 0, for every S > 0, there exists u : [0, T] — > 
[—5,(5] (resp. piecewise constant with value in U n [0,(5]) 
such that for every a in [0, 1] and (j, k) in {1,2} x {1, 2}, 



< e. 



|^,T^ fe )|-|(^-,T^ fe }| 

In other words, up to an arbitrary small error e, it is pos- 
sible to steer the eigenvectors <\>\ and </> 2 , simultaneously for 
every a, to a target with prescribed modulus of coordinates 

on (j>i and (/> 2 . 

The contribution of this note relies on the very same idea 
as [13], namely the computation of finite dimensional Lie 
brackets and a polynomial interpolation. The only novelty 
is that all the steps of the proof come along with explicit 
estimates, which allow to let the dimension of the finite 
dimensional systems tend to infinity and eventually prove 
the infinite dimensional result. 

The main improvements of this work with respect to the 
cited references are 



• the possibly infinite dimension of the ambient space H; 
m the possibility for the spectrum of A to have a contin- 
uous part; 

« the possible (finite or not) degeneracy (or multiplicity) 
of the eigenvalues of A; 

• (in the case of Assumption |2]i the possible unbounded- 
ness of operator B with respect to A, that is, in a case 
where Kato-Rellich theorem does not apply to A + uB. 

E. Content of the paper 

The core of the proof of Proposition Q] is a constructive 
approximate controllability result about the propagator Xf^? 
in some finite dimensional subspaces C n of H proved in 
Section [TT] The precise estimates of Section IH-CI allows 
to let the dimension of C n tend to infinity and eventually 
to prove, in Section IIII-B1 the infinite dimensional result 
for systems satisfying Assumption [2] In Section IIII-C1 we 
will see that the convergence process used for the proof 
of section IIII-BI is actually robust enough with respect to 
perturbation of the spectrum of A to ensure convergence also 
for systems satisfying Assumption [T] The results are applied 
to the example of the 3D rotation of a collection of linear 
molecules in Section HVl 

II. FINITE DIMENSIONAL PRELIMINARIES 

A. Notations and result 

Let N in N, AW^C' be two matrices in u{N) (that 
is, A( N ) T + A iN ) = #) T + = 0). We consider the 



continuum of A^-dimensional systems 

x 1 = A {N) +uaB^x, a G [0, 1] 



(4) 



where x is a point in C N endowed with its canonical Hilbert 
structure (-,-,). For every locally integrable function u, we 
define XT^t the propagator associated with 

We assume that A^ is diagonal in (4>j)j<N, the canoni- 
cal basis of C w , we denote with (— iAj )j<jv the diagonal 
of AW and with b jk := (^,BW^),1 < j,k < N 
the entries of B^ N \ For every j < N, we define TTj, the 
orthogonal projection of C N to Cj = span(0i, . . . , cj>j). 

Proposition 2: Assume that (AW, flW, Ai, (j>x, A 2 , <£ 2 ) 
satisfies Assumption!!] Let f : a G [0, 1] h-> f a £ U(C 2 ) be 
a continuous curve of unitary operators on £ 2 that satisfies 
T° = Idx 2 . Then, for every e > 0, for every S > 0, there 
exists u : [0,T] -> [-6,5] such that, for (j,k) in {1,2} 2 , 



< e. 



|(0 3 ,T^»|-|(^,T^ fc )| 

The proof of Proposition |2] is split in two steps. In a first 
time, after a suitable change of variable, we introduce a 
continuum of two-dimensional auxiliary systems in Section 
IH-BI Classical Lie groups technique, and the associated uni- 
form convergence estimates, to prove approximate ensemble 
controllability of these systems. In a second time, in Section 
IH-CI we use classical averaging techniques to show that the 
trajectories of the systems introduced in Section Hl-B I can be 
tracked, with arbitrary precision, by the system (0}. 



B. An auxiliary system 

We consider the continuum of control systems in U(2) 



a 



hi 
b 2 \e~ 



b 12 e h 

b 22 



a G fO, 11 



(5) 



with initial condition x a (0) = h and control function 9 : 
R — » R. For every piecewise constant function : R — » R, 
for every a, we denote with y e a the propagator of (0. 

Proposition 3: Let f : a G [0, 1] i-> f a G S , C/(£ 2 ) be 
a continuous curve of unitary operators on C 2 that satisfies 
T° = Id£ 2 . Then, for every e > 0, there exists w : [0, T] — > 
[— 7r,7r] piecewise constant such that, for (j,k) in {1,2} 2 , 



K0 3 -,i?' a &>|-K&,T«> fc )| < 

Proof: Let e > 0. There exists a continuous function 
ti : a 4 ii° 6 su(£2) such that, for every a in [0,1], 
|jexp(€> Q ) — T"|| < e. By density of odd polynomials 
mapping, for the norm of uniform convergence, in the set of 
odd continuous functions, there exists a polynomial mapping 
P P a = J2? =0 a 2l+1 Z 2l+1 , with Z u Z 3 , . . . , Z 2N+1 

in su(2) such that ||P Q — v a \\ < e for every a in [0, 1]. 

Lemma 4: Let X and Y two matrices in su(2), and 
(Cj(X, Y))i<j<p a sequence of iterated brackets of X and 
Y. We denote with lj the length of the bracket Cj(X,Y) 
(the length of [X, Y] is 1). Then, for every real sequence 
(/3j)i<j< p , for every T in R, for every e > 0, there exists 
a finite sequence (tk)i<k<m in R such that, for every a in 
[0,1], ||P Q - e TzZU^ aTjC ^ x -y)\\ < e , where P a is the 
product of matrices e tia -V 2Qy • • • e tm - iaX e tmClY . 

Proof: This result is very classical when a = 1 
(i.e., one considers one system only). The uniform version 
presented here (with a in [0, 1]) is basically contained in 
[17]. Because of its importance for our purpose, we give 
below a sketch of the proof of the result. 

We first assume that p = 1 and we proceed by induction 
on the length of C\. From the Baker-Campbell-Hausdorff 
formula, we deduce that, for every 2x2 matrices X, Y, 
there exists a function gx.Y : R — > 0^ tending to 02 X 2 at 
such that, for every t in R, for every a in [0, 1], 



taX taY -taX -taY 

e e e e 



^a 2 t 2 [X.Y]+9x,Y(at) a 2 t 2 



As a consequence, for every 2x2 matrices X and Y, 

= 0, (6) 



lim i 

t-yo t 2 



e taX e taY e -taX & -taY _ e a 2 t 2 [X,Y] 



the convergence being uniform with respect to a in [0, 1]. 
Recall that, for every V, W in su(2), for every n in N, 

\\v n -w n \\ = v{v n - x - w 11 - 1 ) + (V - W)W n ~ 1 

< wv*- 1 - w^w + 1| v- w\\ 

< ra||V-W||. (7) 



Hence, 



I taX taY -totX -taY\ 

e e e e 



^na 2 t 2 [X,Y] 



< 



n e 



taX taY -taX -taY _ a 2 t 2 [X,Y] 



Choosing n = T/t 2 and letting n tend to infinity (and hence 
t tend to zero) gives the result for p = 1, j3i = 1 and l\ = 1. 
The proof for l-y > 1 is very similar, replacing X and Y by 
the suitable iterated brackets in (0. 

A consequence of Zassenhauss formula is that, for every 
2x2 matrices U, V, there exists a locally Lipschitz function 
g : gl(2) x gl(2) x R —> R that vanishes as soon as one of 
its entries vanishes such that, for every t in R, for every a 
in [0, 1], for every j, k in N, 

\\ e tu e tv _ e t(u+v)\\ < t 2 g{UjV ,t). 

The proof of Lemma |4] for p > 1 and j3j not necessarily 
equal to 1, follows by choosing t = T/n for n large enough 
and using once again (O. ■ 
We come back to the proof of Proposition |3] After the 
time dependent change of variable 



exp 



—tot 



foil 






b 22 



(9) 



the system (0 reads 



b 12 e" 
fe 2 ie~ ie+ * (bll " b22) 



-t(bi 



or y' a = M£y a , defining u = 0- it(b n - b 22 ), with 



Ml = a 







bi 2 e w 



b 21 e- w 



(10) 



(ID 



Lemma 5: For every </> in C 2 , for every t in R, for every 
locally integrable 9 : R — > R, the moduli of the coordinates 
in the canonical basis (</>i, (j> 2 ) of C 2 of y a (t)ip and x a (t)ip 
are the same. 

Proof: From (O, the coordinates of y a and x a are 
equal, up to a phase shift depending on time and a. ■ 
Thanks to Lemma 2J Proposition [3] follows if, for every /, 
the matrix Z 2 i+i defined above can be realized as a linear 
combination (with real coeeficients) of brackets of length 
exactly equal to 21 + 1 of the matrices M^,v G [— vr,7r]. 
Notice that 



and M 2 



1 



i&i2 



Straightforward computations give, for every k in N, 

■21 „ ) • ( 12 ) 



Ml a 



&i 2 |6i 2 | 2fe 



6 2 i|foi 2 | 2fc 

Proposition [3] follows from the fact that b\ 2 ^ 0. ■ 

C. Averaging techniques 

We define the AT x N matrix by N®(j,k) = 
for every j, k in {1,... ,A} 2 but A^(l,2) = ab 12 e w and 
Af (2, , 1) = -A»(l, 2). In particular, N e a belongs to su(A) 
and 7r 2 A^7T2 = M e a . 

Let us come back to the proof of Proposition [2] From 
Proposition |3] it is enough to show that, for every 9,t in R 
and every e > 0, there exists u £ : [0,T e ] — > (— <5, (5) such 
(8) that, for every a in [0, 1], W^X^T,) - e tN *\\ < e. This 



is exactly the content of Proposition [6] whose proof is given 
in [20]. 

Proposition 6: Let u* : R + — > R be a locally integrable 
function. 



Assume that u* is periodic with period T = 



■2 IT 



A2-A1I 



and 



that J u*(r)e i(Al Am)r dr = for every {l,m} such that 

{l,m}°n {1,2} ^ and A; - A m G (Z \ {±l})(Ai - 
A2) and bi, n 7^ 0. For every n, define v n : t 1— > 
1/n J Q |w*(s)|ds and the TV x TV matrix with entries 

mj,k = V/o rM *( s ) ei(A2 " Al)Sds //o T |"*( s )l ds - 

If / M*(r)e l(A2_Al)T dr ^ 0, then, for every n in N, for 
Jo 

every i < nT , 



J(C+1)||BW|| 



(13) 



with 



2|6i, 2 | / V(r)e i (^i-^)T dr 



J= / |u*(r)|dr, 
/o 



K 



IT* 



,C 



sup 

0,fc)€A 



sm 



|Aj-A fc | 
A2 — Ai 



where A is the set of all pairs (j, k) in {1, . . . , A} 2 such that 
b jk ^ and {j, k} n {1, 2} ^ and |A, - A fe | £ Z|A 2 - \ x \. 

Proof: (Proposition [2]i We apply Proposition [6] with 
w* periodic with period T = \\^-x^\ an( l satisfy- 
ing f M*(r)e i(A2 - Al)r dr ^0 / u*(r)e i(A ^ Am)T dT = 

Jo Jo 
for every {Z, m} such that {I, m}n{l, 2} ^ and A; — A m € 

(Z \ {±l})(Ai — A2) and bi m ^ 0. Such a it* can be chosen 

of the form t i-> cos((A2 — Ai)t — 9) or piecewise constant 

with value in {0, 1} (for an explicit construction of such a 

function, see [12]). 

To ensure that e v ™ XI (*) Mt is e-close to e tM °, defined as 

in ( fTTT ). one chooses t such that v n (t) = b\2& lB /m\ 2 < nT*. 

One can check from the definition of v n that t/(nr*) tends 

to 1 as n tends to infinity, where r* is defined by 

r* = ( *€-\ . 



'12 



The final step in the proof of Proposition [2] is to get rid 
of the phase e tA in estimate dO} . We use the Poincare 
Recurrence Theorem with the mapping 



R 



R n+1 /Z™ +1 



s y y [ tt^s. -7^- s . 

\ 27T ' 2"7T > 



2-7r ' r* / 



on the rt + 1 dimensional torus. For every e > 0, there exists 
a sequence (sk)keisi that tends to infinity such that s^Aj is 
e close to 27rZ and sj. is e-close to r*Z. The sequence of 
controls \r* / Sk\u* gives Proposition [2] by letting k tend to 
infinity. ■ 



III. INFINITE DIMENSIONAL ESTIMATES 

A. Heuristic of the proof 

In this Section, we proceed to the proof of Proposition 
[U Inspired by Section [II] it is enough to show that the 
projections of each of the infinite dimensional systems (O 
can track, with an arbitrary precision, the trajectories of the 
2x2 system ©. 

To begin with, we consider in Section IIII-BI a system that 
satisfies Assumption [2] The proof is a uniform version of the 
Section 4 of [20] which is valid for one particular a. 

To prove Proposition[T]for systems that satisfy Assumption 
[U we first estimate the robustness of the results of Section 
IIII-BI against a perturbation of the spectrum of A. The 
conclusion will follow from the Von Neumann approximation 
theorem. As in Section IIII-BI the method of the proof 
in Section IIII-CI is similar to the one used in [20], the 
only difference lying once again in the uniformity of the 
convergence estimates with respect to a in [0, 1]. 

B. If the eigenvectors of A span a dense subspace of H 

Let (A, B, U, A, $) satisfy Assumption^ 9 in [— tv, it] and 
r, e > 0. We aim to find a periodic control u* with period 
2ir/ (A2 — Ai), n in N and T in R + such that, for every a 
in [0,1], 



T 



u* /n,a 
T,0 



< e. 



Since <fii and </> 2 belong to the domain of B, the sequences 
(&i,j)i£N and (&2,z)igN are in £ 2 . Hence, there exists N in 
N such that \\tt 2 B(1 - n N )\\ = ||(1 - 7nv)-B7r 2 || < 5e/(2r). 
Define ui = A 2 — Ai and u* with period 2it/uj in such a 
way that the N x N matrix M> of Proposition [6] is equal 
to M{ and the efficiency | ffi u*(s)e iujs ds\/ \u*(s)\ds 
of u* for the transition (1,2) is larger than 2/5. This can 
be done, for instance, with t M> cos(wt — 9) (efficiency 
7r/4) in the case where B is bounded or, in the general case 
of Assumption [2] with a piecewise constant function taking 
value in {0, 1} as described in [12]. 

For a given n to be precised later, consider system (|3} 
with control u n = u* Jn in projection on span(</>i, . . . , (J)n)' 

TTN^r^h = (AW +u n {t)aB^)T u r a ^ 

+ U n (t)TT N aB(l - 7TAr)T^- Q ^.(14) 

From the variation of the constant, we get, for j = 1,2, 



+/ u n (s)A ( ^ (t, s)n N aB(l - n N )T t n ' (f>jds.(15) 



Project ([T5l l on span(0i, fa), and recall that ttn^ 
TT2TTN — 7i"2 for N > 2: 



(16) 



+ j u n (s)n 2 X^ (t, s)ir N aB{l - n N )T t ' 1 ' fads. 



Define, for every t, s in R, the bounded linear mapping 

[lT2,Xffif(t,s)] := 7T 2 O X"™!™(t,s) - X^ a (t, S ) O 7T 2 . 

Equation ( TToT ) reads, for j = 1,2, 

Tt 2 r t n ' a tpj - TT2X^ a {t,Q)<f)j = 



u n (s)X^ a (t, s)ir 2 aB(l - 7r Ar )T""' Q ( ?!) i ds (17) 



+ J u n (s)[ir 2 , X^™'" (t, s)]Tr N aB(l - 7r A r)T"" ,Q (/) i ds. 

Extend the definition of Mg to H by Aft = on Cj^ and de- 
fine the linear operator E% a (t) := X ( %(i,0) - e"» 1] W M »! 
Since the commutator [tt2, M@] = -kzMq — MS 1 ^ vanishes, 
we have, for every t in R, 



\\[7r 2 ,X^(t,0)}\\ = \\[«2,e« [ - 1] ®^+E{$(t)]\\ 
= ||[7r a ,£^(*)]||<2||^(t)||. 

Note also that, for every t in R, 

|| [Tra, X ( u ^(0, t)]\\=\\X$fp, t)[X$fi, 0), TrajX^O, t)\\ 
< 2||£^(t)||. 

For every s, t in R, 

[7r 2 ,X ( ^ } (t, S )] 

= K 2 Xfc (t, 0)Xfc (0,8)- Xfc (t, 0)Xfc (0, S )tT 2 

= X^ftfi) [n 2 ,X^ (0, a)] + [tt^X^ (t, 0)}X^ } (0,s). 
Finally, we get, for every (s, t) in R 2 , for every n, N in N. 

<4||^(t)||. (is) 

From ((17) and (fTHT ). since ||7r 2 B(l - 7Tjv)|| < e/K, 



|k 2 T-(<)^ 2 -7r 2 X ( ^(t,0V 2 || 

< £ + 4||£^(t)||Ar||7r iV -B(l-7r J v)||. (19) 

From ( fT3l >. sWt<v n {K) ll-®o\n(*)ll tends to zero as n tends 
to infinity. For n large enough, for every a in [0,1] and 
* < v n (r), 

H^W^H - 4^11^5(1-^)11- 

Proposition [6] completes the proof of Proposition Q] in the 
case where (A, B, U, A, <£>) satisfies Assumption [2] 

C. If A has a mixed spectrum 

Assume that (A, B, Ai, (/)%, A 2 , ^2) satisfy Assumption [1] 
From Theorem 2.1, page 525, of [21], for every 77 > 0, 
there exists a skew-adjoint operator A v such that admits 
a complete family of eigenvectors (& v ) associated with the 
family of eigenvalues (A v ), A<fi = A v cj> for every eigenvector 
4> of A and ||^4 - A,,| < 77. 

For every locally integrable u, we denote with T™ the 
propagator of ^i/j = (A n + uB M )ip- 

The scheme of the proof is a follows: the result is known 
(from Section IIII-Bb for the system (A v ,B,ll,A v ,$ 7l ), 
which satisfies Assumption[2] we chose u* : t >-> cos(o;t — 9) 



(this function is the "shape" of the control pulses, it does 
not depend on 77 nor e). For every 77, e > 0, 9 in [— %, tt] 
and r > 0, there exists an integer n n and a positive real T v 
satisfying, for every a in [0, 1], 



T„,0 



_rAf! 1 



< e. 



Notice that, for every £ in R, 1 1 T 



u* J n,ct,r) 
t,0 



' 1 t,0 1 '- 



The crucial point in the proof of Proposition [T] for systems 
satisfying Assumption [T] is the existence of a uniform bound 
on T n , that depends only on r and s, and not on 77. This 
follows from ( fTST l, where the only variable depending on 77 
is 



C = sup 
(i,fc)eA 



J T u*(r)e i ^- Ai! ) r d7 



sin (tt^ 



-A*, I 



where A is the set of all pairs (j, k) in {1, . . . , A^} 2 such that 
b jk and {j, k} n {1, 2} ^ and \X 3 - X k \ £ Z|A 2 - Ai|. 
(Notice that ||B (JV) || is bounded, for every N by ||B||.) 

Straightforward computation gives C < 2/d where d 
is the distance of the set {2Ai — A 2 ,Ai,A 2 ,2A 2 — Ai} to 
the continuous part of the spectrum of 1A. This distance is 
not zero by Assumption 11 161 what concludes the proof of 
Proposition Q] 

IV. EXAMPLE: ROTATION OF A MOLECULE 
A. Modeling 

The description of the physical system we consider is a 
toy model inspired by the physical system described in [22]. 
It has already been thoroughly studied, see for instance [23], 
[12] or [24]). We consider a polar linear molecule in its 
ground vibronic state subject to a nonresonant (with respect 
to the vibronic frequencies) linearly polarized laser field. The 
control is given by the electric field E = u(t)(E\, E 2 , E3) 
depending on time and constant in space. We neglect in this 
model the polarizability tensor term which corresponds to 
the field-induced dipole moment (see for instance [25] or 
[26]). 

Let V be a fixed plane in the space. We are interested 
in the orientation of the orthogonal projection of a set of 
molecules in the plane V (given by one angle, in contrary to 
the orientation of the molecule in the space which is given by 
two angles). We neglect the interaction between molecules, 
and consider only the interaction between the molecules and 
the external field. Our aim is to control the orientation of 
projection of the molecule in V, whatever the angle of the 
molecule could be with V. 

Up to normalization of physical constants (in particular, 
in units such that h = 1), the dynamics of each molecule is 
ruled by the equation 



. dfl>(9,t) 
1 dt 



-Aip + ui(t) cos 9 sin ipip (9, (p,t) (20) 



where 9 is the angular coordinate in V and ip is the angle of 
the molecule with V ', which is assumed to be constant for 
the sake of simplicity, A is the Laplace-Beltrami operator 
on the circle S = R/27rZ, The wavefunction ip(-,t) evolves 



in the unit sphere S of H = L 2 (S, C) endowed with scalar 
product (f,g) = f(s)g(s)ds. 

The operator A = iA is skew-adjoint in H, with domain 
H 2 (S,C) and has discrete spectrum. Define <fio : 9 i— > 
l/-\/27r and, for every fc in N, <p2k-i '■ & l— ^ cos(fc0)/y / 7r 
and 02fc : n- sm(k0) / y/n . The two functions <p2k+i and 
(/>2fe are eigenvectors of A, associated with eigenvalue — ifc 2 . 

The operator B : ip M> — icos(6)ip is bounded. Straight- 
forward computations show that \{(f>o, B, = l/v2 and 
(<f>j,(f>k} — if the parities of j and fc are different or if 
\j -k\>2. 

B. Result 

Assume that a bunch of molecules is in the state (j>o at 
t = 0. We aim to transfer to the state tj>\ all the molecules 
for which ip > tt/3 and to keep all the molecules for which 
ip < 7r/6 in the state 4>q. 

From Proposition Q] applied to (A, B, 0, <f>o, 1, ^2) which 
satisfies Assumption Q] this is possible, up to the phase and 
to an arbitrary small error e > 0. 

V. CONCLUSIONS 

A. Some comments on the result 

While our construction is completely explicit (simple 
formulas are available for the control laws and come along 
with precision and time estimates), the convergence toward 
the target is extremely slow and cannot be used for actual 
control of real systems. This well-known fact is due to the 
very poor efficiency of tracking strategies via Lie brackets. 

B. Perspectives 

The presented results may certainly be improved in many 
ways. For instance, the author conjectures that it is possible 
to replace T in Proposition \T\ by a unitary transformation 
of Cn with N > 2 or to extend the result to systems for 
which the free Hamiltonian A has a mixed spectrum and the 
coupling Hamiltonian B is unbounded. 

VI. ACKNOWLEDGMENTS 

This work has been partially supported by INRIA Nancy- 
Grand Est, by French Agence National de la Recherche ANR 
"GCM" program "BLANC-CSD", contract number NT09- 
504590 and by European Research Council ERC StG 2009 
"GeCoMethods", contract number 239748. 

References 

[1] J. M. Ball, J. E. Marsden, and M. Slemrod, "Controllability for 

distributed bilinear systems," SIAM J. Control Optim., vol. 20, no. 4, 

pp. 575-597, 1982. 
[2] G. Turinici, "On the controllability of bilinear quantum systems," in 

Mathematical models and methods for ab initio Quantum Chemistry, 

ser. Lecture Notes in Chemistry, M. Defranceschi and C. Le Bris, Eds., 

vol. 74. Springer, 2000. 
[3] K. Beauchard, "Local controllability of a 1-D Schrodinger equation," 

/. Math. Pures AppL, vol. 84, no. 7, pp. 851-956, 2005. 
[4] K. Beauchard and J.-M. Coron, "Controllability of a quantum particle 

in a moving potential well," J. Fund. Anal, vol. 232, no. 2, pp. 328- 

389, 2006. 



[5] V. Nersesyan, "Global approximate controllability for Schrodinger 

equation in higher Sobolev norms and applications," Ann. Inst. H. 

Poincare Anal. Non Lineaire, vol. 27, no. 3, pp. 901-915, 2010. 
[6] , "Growth of Sobolev norms and controllability of the 

Schrodinger equation," Comm. Math. Phys., vol. 290, no. 1, pp. 371- 

387, 2009. 

[7] K. Beauchard and V. Nersesyan, "Semi-global weak stabilization of 
bilinear Schrodinger equations," C. R. Math. Acad. Sci. Paris, vol. 348, 
no. 19-20, pp. 1073-1078, 2010. 

[8] K. Beauchard, J. M. Coron, M. Mirrahimi, and P. Rouchon, "Implicit 
Lyapunov control of finite dimensional Schrodinger equations," Sys- 
tems Control Lett, vol. 56, no. 5, pp. 388-395, 2007. 

[9] M. Mirrahimi, P. Rouchon, and G. Turinici, "Lyapunov control of 
bilinear Schrodinger equations," Automatica J. IFAC, vol. 41, no. 11, 
pp. 1987-1994, 2005. 
[10] M. Mirrahimi, "Lyapunov control of a particle in a finite quantum 
potential well," in Proceedings of the 45th IEEE Conference on 
Decision and Control, December 2006. 
[11] T. Chambrion, P. Mason, M. Sigalotti, and U. Boscain, "Controllability 
of the discrete-spectrum Schrodinger equation driven by an external 
field," Ann. Inst. H. Poincare Anal. Non Lineaire, vol. 26, no. 1 , pp. 
329-349, 2009. 

[12] U. Boscain, M. Caponigro, T. Chambrion, and M. Sigalotti, "A weak 

spectral condition for the controllability of the bilinear Schrodinger 

equation with application to the control of a rotating planar molecule," 

Comm. Math. Phys., vol. 311, no. 2, pp. 423^55, 2012. 
[13] J.-S. Li and N. Khaneja, "Control of inhomogeneous quantum 

ensembles, " Phys. Rev. A, vol. 73, p. 030302, Mar 2006. [Online]. 

Available: http://link.aps.org/doi/10. 1 103/PhysRevA.73.030302 
[14] K. Beauchard, P. S. Pereira da Silva, and P. Rouchon, 

"Stabilization for an ensemble of half-spin systems," Automatica 

J. IFAC, vol. 48, no. 1, pp. 68-76, 2012. [Online], Available: 

http://dx.doi.Org/10.1016/j.automatica.2011.09.050 
[15] K. Beauchard, J.-M. Coron, and P. Rouchon, "Controllability issues for 

continuous-spectrum systems and ensemble controllability of Bloch 

equations," Comm. Math. Phys., vol. 296, no. 2, pp. 525-557, 2010. 

[Online]. Available: http://dx.doi.org/10. 1007/s00220-010- 1008-9 
[16] J.-S. Li and N. Khaneja, "Ensemble controllability of the bloch 

equations," in Decision and Control, 2006 45th IEEE Conference on, 

dec. 2006, pp. 2483 -2487. 
[17] , "Ensemble control of Bloch equations," IEEE Trans. Automat. 

Control, vol. 54, no. 3, pp. 528-536, 20 09. [Online]. Available: 

http://dx.doi.org/10.1109/TAC.2009.2012983 
[18] Z. Leghtas, A. Sarlette, and P. Rouchon, "Adiabatic passage and 

ensemble control of quantum systems," Journal of Physics B, vol. 44, 

p. 154017, 2011. 

[19] A. Zlotnik and S. Li, "Iterative ensemble control synthesis for bilinear 

systems," in Decision and Control (CDC), 2012 IEEE 51st Annual 

Conference on, dec. 2012, pp. 3484 -3489. 
[20] T. Chambrion, "Periodic excitations of bilinear quantum systems," 

Automatica J. IFAC, vol. 48, no. 9, pp. 2040-2046, 2012. 
[21] T. Kato, Perturbation theory for linear operators, ser. Classics in 

Mathematics. Berlin: Springer- Verlag, 1995, reprint of the 1980 

edition. 

[22] U. Boscain, M. Caponigro, and M. Sigalotti, "Controllability of the 
bilinear schrodinger equation with several controls and application to a 
3d molecule," in Decision and Control (CDC), 2012 IEEE 51st Annual 
Conference on, dec. 2012, pp. 3038 -3043. 

[23] U. Boscain, T. Chambrion, P. Mason, M. Sigalotti, and D. Sugny, 
"Controllability of the rotation of a quantum planar molecule," in 
Proceedings of the 48th IEEE Conference on Decision and Control, 
2009, pp. 369-374. 

[24] N. Boussai'd, M. Caponigro, and T. Chambrion, "Periodic control laws 
for bilinear quantum systems with discrete spectrum," in Proceedings 
of the American Control Conference, 2012. 

[25] M. Morancey, "Explicit approximate controllability of the schrdinger 
equation with a polarizability term," Mathematics of Control, 
Signals, and Systems, pp. 1-26, 2012. [Online]. Available: 
http://dx.doi.org/10.1007/s00498-012-0102-2 

[26] N. Boussaid, M. Caponigro, and T. Chambrion, "Approximate con- 
trollability of the Schrodinger equation with a polarizability term," in 
Proceedings of the 51st IEEE Conference on Decision and Control, 
december 2012, pp. 3024 -3029. 



